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Q_i' Abstract 



We compute the two-loop electroweak corrections to the radiative decays of the B 
meson in the SM. Electroweak effects reduce the Wilson coefficient C% (M w ) by 2.6% 
for a light Higgs boson of about 100 GeV and are less important for a heavier Higgs. 
The leading term of a heavy top expansion of our result differs from the one obtained 
in the gaugeless approximation where only top quark Yukawa couplings are considered: 
we discuss the origin of the discrepancy and provide a criterion for the validity of the 
gaugeless approximation. As a byproduct of the calculation we also obtain the O(a) 
corrections to the Wilson coefficient of the four-fermion operator Q 2 - A careful analysis 
of the interplay between electroweak and QCD effects leads to an overall 2% reduction of 
the total branching ratio for B — > Xg'y due to purely electroweak corrections. For a light 
Higgs boson, the up-to-date SM prediction is £>i? 7 = 3.29 x 10~ 4 . 



1. Radiative B decays represent one of the most important probes of new physics and a 
major testing ground for the Standard Model (SM). They already place severe constraints 
on many new physics scenarios. The present experimental accuracy for the branching ratio 
of B — > X s 7 (BR 7 in the following) is about 15% and is expected to improve significantly 
in the near future, both at CLEO and at the B factories. 

On the theoretical side, since precise predictions in the SM are particularly important, 
the subject has reached a high degree of technical sophistication. Indeed, perturbative 
QCD corrections are very sizeable |2| and give the dominant contribution; they are best 
implemented in the framework of an effective theory obtained by integrating out the heavy 
degrees of freedom characterized by a mass scale M > M w . At lowest order in this 
approach the FCNC processes B — > X s 7 and B — > X s g proceed through helicity violating 
amplitudes induced by the magnetic operators 

Q 7 = £^rn b s L a^b R , Q 8 = ^m b s L a^t a b R G% . (1) 

A few years ago the renormalization group improved QCD calculation has been completed 
at the next-to-leading order (NLO) [§-0, reducing the uncertainty from uncalculated QCD 
higher orders to about 5%. More recently, NLO predictions have been made available in 
some new physics models as well [0.0]. 

There has also been progress concerning QED and electroweak radiative corrections: 
after Czarnecki and Marciano considered all the leading QED logarithms || , their interplay 
with QCD corrections has been studied in H[TI|, an d Strumia JTT| has calculated the 



leading term of the Heavy Top Expansion (HTE) of the electroweak two-loop corrections 
using the gaugeless limit of the SM. As for non-perturbative effects, they seem to be under 
control [0 , although some aspects may still need a more detailed investigation |13| . 



In addition to uncalculated radiative corrections and non-perturbative effects, the cal- 
culation of BR 7 is also affected by the uncertainties on the input parameters (the CKM 
matrix elements, the semileptonic branching ratio BR^^, etc.). In fact, the latter bring the 
overall theoretical error to about 10%. As the parametric uncertainties (especially those on 
the CKM elements, a s , BRsl, and M t ) are expected to decrease soon, and given the crucial 
importance of this decay mode, it seems appropriate to try and refine the SM prediction 
as much as possible. 

In this note we reexamine the two-loop electroweak contributions to radiative B decays 
and present the result of a calculation where only some photonic effects have been neglected. 
Moreover, we update the SM prediction of BR 7 using the latest experimental inputs. The 
final result is expressed by a compact formula that summarizes the dependence on the 
input parameters. 

2. Although generally small, two-loop purely electroweak effects are sometimes very 
important: an example is provided by the precision observables of the SM, like the ef- 
fective sine measured on the Z° pole and the W mass, where radiative corrections up 
to 0(g 4 Mt /M^) |14|,[15|] are now routinely included in the analysis with important con- 



2 



sequences in the electroweak fits fT5, 16||. Moreover, by fixing the normalization of the 
electroweak coupling, two-loop effects reduce the electroweak scheme dependence of the 
SM prediction, which can be quite large — also for FCNC processes [17] . 



As mentioned above, in the calculation of BR 7 the leading large logarithms of QED 
origin are now under control, as a resummation of all aa™" 1 (In rrib/M w ) n terms has been 
completed USHlOfl. Apart from that, our knowledge of electroweak effects in b — > s'j is 



limited to the subset of two-loop fermion loop corrections calculated in [g] and to the 
leading term of the HTE of pl| . In fact, the two results are numerically very different 
- about —2.3% and less than —0.7%, respectively, on the Wilson coefficient at M w . The 
leading term of the HTE was calculated in []TTJ] using the gaugeless limit of the SM, i.e. in a 



Yukawa theory where the heavy top couples only to the Higgs doublet, setting M w = and 
keeping the Higgs mass M H finite and arbitrary. In the presence of external gauge bosons, 
these can be considered as background sources. This approach presents a few limitations 
that also motivate our new calculation: 

• the lowest order contribution to the Wilson coefficient of is a function of the top 
mass whose HTE converges very slowly. Using x t = M^/M^ ^4.7 and writing 
explicitly the numerical values of the successive 0(l/x™) terms, it reads 

(o), , _ x t (7 - 5x t - 8xf) x 2 t (3x t -2)\nx t _ 
° 7 {Xt) ~ 24(x,-l)3 + 4(0,-1)4 " (2) 

= -- - 0.010 + 0.070 + 0.046 + 0.021 + ... = -0.195 
3 

where the ellipses represent contributions 0(1/ x^) or higher. The leading HTE is 
therefore unlikely to provide anything more than an order of magnitude estimate of 
the two-loop electroweak contribution. In this respect, the similar case of B — B 
mixing |I7| is very instructive: for realistic values of the top mass the complete two- 
loop electroweak correction is not well approximated even by the first three terms of 
the HTE and the leading HTE term is numerically far from the complete result. 

• even assuming the leading HTE term to be representative, it should not be expected 
to give an accurate result for a light Higgs mass, M H « 0(M W ), because it is obtained 
by setting M w = fl5|]. On the other hand, present electroweak fits show a decisive 



preference for a light Higgs boson, M H < 215 GeV at 95% C.L. fllSf. 



the gaugeless limit has often been used to compute the leading HTE term, but it is 
known (T^j that in some cases it does not reproduce the correct result. In the following 



we explain why it fails for radiative B decays and provide a general criterion for its 



use. 



A complete calculation of all electroweak effects in radiative B decays in the frame- 
work of effective Hamiltonians is a very complex enterprise which involves other operators 
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in addition to those of Eq. (Tl]). In fact, the analysis should be aimed at resumming all 
«a" ()iamb/M w ) n effects. The procedure is summarized, for instance, in Hl9| . Its necessary 
steps would be: (i) the calculation of two-loop 0(a) matching conditions for Q7,8 at some 
0(M W ) scale — this involves also their QED mixing with all other operators — and of the 
0(a) contributions to various four quark operators; (ii) QED-QCD running of the Wilson 
coefficients to the B mass scale — this would require a three loop computation of the 
anomalous dimension matrix similar to that of ||; (iii) calculation of the one-loop QED 
matrix elements of the various operators — the determination of these matrix elements 
depends sensitively on the precise experimental conditions. 

An important simplification can be obtained by keeping only the first term in an ex- 
pansion around s w = sin^ = 0. This is equivalent to considering a SU(2)l theory with 
a background photon field and removes all the light virtual degrees of freedom. In partic- 
ular, all the diagrams with virtual photons and all infrared (IR) divergences drop out of 
the two-loop calculation in a gauge-invariant way. Step (i) is therefore much simpler as 
the calculation of the two-loop b — > 37 and b — > sg amplitudes gives us directly the scheme 
independent 0(g 2 ) correction to C 7>8 , respectively (g = e/s w is the SU(2) L coupling). 
Moreover, this simplification avoids completely steps (ii) and (iii), because they are both 
driven by purely photonic effects suppressed at least by Q u \Qd\ s w ~ 0.05 with respect to 
pure SU(2) L contributions. 



In analogy to [19], we complement this approximation scheme by keeping also the 
0(g 2 M 2 /M^) contributions that vanish as s w — > 0. In practice, we therefore expand the 
two-loop b — > s r y(g) amplitude in powers of s^: 

'A + A lS 2 w + O(s 4 w )} (3) 



A^ = g* 



and retain only Aq and the 0(M 2 ) part of A%. This is likely to be a sufficiently good 
approximation, as suggested by those cases [|rF|>E9 where it has been possible to compare 
it with the complete SM result. We recall that Ap |2(J and Bq — Bq mixing involve 



amplitudes conceptually similar to those under consideration: forbidden at tree level and 
induced at one-loop by virtual SU(2)l effects. Later on, we will give an estimate of the 
residual uncertainty. 

We calculate analytically the two-loop amplitudes in the Feynman background gauge 
with a background photon (gluon). A few thousand diagrams are automatically generated 
by the package FeynArts 2.2 (the topologies are shown in Fig. 1). After setting to zero 
all light fermion masses but the b quark mass, they can be reduced to a few hundred equiv- 
alence classes, which we have actually computed. Due to the GIM mechanism, the CKM 
cofactor of each equivalence class is always proportional to A t = V t bV t * s . The extraction 
of the magnetic penguin amplitude and the two-loop integration are performed as in H . 
All the steps of the calculation have been implemented in two independent and completely 
automatic codes that involve various combinations of Mathemathica [^] and Form 
routines. Although the result can be expressed in terms of logarithms and dilogarithms, it 
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Figure 1: Two-loop topologies for b 



S7. 



is rather lengthy and we will present instead accurate numerical approximations. 

A peculiarity of the two-loop calculation for these processes is the presence of diagrams 
containing anomalous fermionic loops (triangles). It is well known that the naive definition 
of anticommuting 75 in n dimensions that we employ in the rest of the calculation fails for 
these diagrams because it leads to algebraical ambiguities and cannot reproduce the axial 
anomaly. Our solution consists in calculating Dirac structures containing an odd number 
of 75's — i.e. those leading to the anomalous term — using anticommuting 75 in strictly 
four dimensions, which is possible because of their apparent UV convergence. The anomaly 
cancellation then guarantees the absence of both anomalous and ambiguous terms in the 



sum of all diagrams (see [p5| , app. C). We have also checked our results for these specific 
terms using the HV definition of 75 in n dimensions [ 24f . From a formal point of view, the 



equivalence of the two methods follows from the absence of non-invariant counterterms for 
the odd 75 part in the HV case [ 25| , as can be also seen in full generality using the powerful 



formalism of 26 



The renormalization is performed following the simple framework of |T7]], to which we 
refer for a detailed discussion (see also |15|) and for the notation. We recall that the top 
mass is renormalized on-shell as far as electroweak effects are concerned, but it is customary 
to use an MS definition for the QCD effects. Although the choice of the scale fit for the 
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MS top mass is a matter of convention, the NLO QCD corrections depend sensitively on 
\i t . For simplicity, we follow [|| and in the numerics we set \i t = \i w = M w and employ 
M t = M t {M w ) = 175.5 ± 5.1 GeV obtained from the pole mass value Mf ole = 174.3 ± 5.1 
GeV . No renormalization of the electric charge is necessary in our approximation. The 
renormalization of the b quark mass, not needed in [ITT 
counterterm reads 



is also performed on-shell. The 



167T 2 \Mt 



3 < 

8e V 



r 



3 + 8x t - hx\ 3 (1 - 3x t + xf) In x t 



16 (x t - 1) 



fa - iy 



+ 0(si 



One should keep in mind that the nib factors in Eq. ([[[) originate either from the 6-quark 
Yukawa coupling or from the use of on-shell equations of motion. In the latter case, nib 
should not be renormalized as it is on-shell by definition. This is the mass appearing 
explicitly in the projector of Eq.(14) of ||. Indeed, besides the magnetic operators of 
Eq. (P, there are additional off-shell operators that project onto Q-j^ when the external 
momenta are set on-shell, i.e. Q w ~ e s L {p, a^ u }bL and the analogous one with gluon 
fields ||. In correcting the external fields, one should take into account that the chirality 
of the h quark is different in Q 7 and Qiq. The external leg corrections of |T7| correspond to 
the correct LSZ factors and implement the renormalization of the CKM matrix according 
to [^] within our approximations. We recall that this gauge invariant definition of the 
CKM matrix is the most appropriate to the present low-energy measurements because, 
unlike an MS renormalization, it avoids 0(g 2 ) corrections not suppressed by GIM and 
proportional to (mf +m|)/(mf — where niij are light quark masses [0. Notice also 
that in our framework 6(Z^)ij = 0. 

According to the standard procedure, we will normalize BR 7 to the semileptonic branch- 
ing ratio, BRsx. This fixes the normalization of the electroweak coupling but requires the 
inclusion of the one-loop electroweak corrections to BR^^ ||28||. It is straightforward to see 



that, up to O(s^) terms that we neglect, these are the same that enter the muon decay. 
Hence, in this respect the use of BRg^ is effectively equivalent to that of the Fermi constant 
measured in muon decays, G M , and the coupling renormalization proceeds as described for 



this case in ||17|| . We incorporate the complete one-loop correction to the muon decay 
amplitude, without taking the s w — > limit, on the ground that this is an independent 
process for which the complete correction is available. Notice that the leading part of the 
photonic corrections to BR^^, characterized by large logarithms and not considered in our 
calculation is part of the 0(aa™ _1 ln n nib/M w ) analysis of |8}-|lT)|] and is included in our 
numerical results. 

We now recall that the regularization scheme-independent quantity entering the calcu- 
lation of BR 7 is not C-r(nb) but a combination Cj^(fib) of this Wilson coefficient and of 
the coefficients of the four fermion operators with mixed chirality [§,^9|. The case of QED 
effects has been also considered in ||10|| . It turns out that, as far as electroweak corrections 
are concerned, the two scheme-independent quantities relevant for B — > X s 7 and B — > X s g 
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where C^g are the Wilson coefficient of two electroweak penguin operators 
| (sd)y-A ^ e q {qq) v+ Ai Q% = | 
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aJV+A- 



(4) 



(5) 



q=u,d,s,c,t 



q=u,d,s,c,b 



Notice that in QCD there are other contributions to Eq. (f|) As C 8 p (M vy ) = at leading 
order and Gf(M w ) is proportional to s 2 ^, we need to consider only the part of Cf(M w ) 
enhanced by M 2 (which actually approximates the full Wilson coefficient very well ||19|| ). 

Our results for the 0(g 4 ) contributions to B — > X s j(g) can therefore be written as 
additional contributions to C^J {M w ) and are accurately approximated by 
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(6) 



where the SU(2) L coupling g can be calculated from the relation g 2 = ^\p2G il M^ v . We 

0.23145 for the 0{g 2 s 2 w M 2 ) contributions that we retain. 



80.419 GeV and s 2 w 



use M w 

Eqs. (|B|) reproduce accurately (within 1%) the analytic results in the ranges 100 < M H < 
250 GeV and 165 < M t < 180 GeV. We stress that Eqs. (|]) are independent of the choice of 
the scale \x t in the QCD top mass definition: it is sufficient to calculate M t [}it) and employ 
it in Eqs. (|6]). Differences between different choices are present in the QCD corrections to 
BR 7 but are higher order effects as far as the present calculation is concerned. 

The numerical relevance of our corrections to C%/(M W ) is shown in Fig. 2 for M t = 
175.5 GeV: at M H = 100 GeV the Wilson coefficients of £^7,8 are reduced, respectively by 
2.6% and 0.7%. As a measure of the uncertainty due to the expansion around s w = we 
use the difference between the complete correction to the muon decay and its s w — > limit, 
which amounts to about 0.5%. This seems to us a realistic estimate of the error due to our 
approximation and we will use it in the following. If we consider only fermionic loops we 
reproduce the results of || for CV, which lead to a —2.3% reduction of Cj(M w ). Although 
purely accidental, the closeness of this fermion loop approximation to our complete result 
for a light Higgs, M H « 100 GeV is impressive. 

3. Let us now consider the HTE of our results and see how it compares with existing 



analyses. In units # 2 /(167r 2 ) M t /(2M, 
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Figure 2: Electroweak corrections to the Wilson coefficients C 7t8 (M w ). The solid lines represent 
our results with their error estimates, the dashed lines their leading HTE, and the dotted lines 
the results of the gaugeless approximation. 
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where h t = M%fM t , Li 2 (x) is the dilogarithmic function, and <p(x) is given in Eq. (48) 



of pT7|| . As can be seen in Fig. 2, the leading HTE term approximates our full result 
very poorly, especially for a light Higgs. We have also studied the convergence of the HTE, 
calculating its first three terms, and found that for realistic M t values they do not converge, 
in a way very similar to |17| . Our Eq. (0) differs from the analogous one in |ll| by a term 



I ( \ — ~3~ ) ~f~' ^ n ^ e °t ner hand, we agree with |TT| if we perform the calculation 



in the gaugeless limit. This is not surprising because it is known [18[ that the gaugeless 
approximation does not always include all leading M 2 contributions. 

To understand better this point, notice that for a asymptotically heavy top both the 
top Yukawa coupling, g t = g M t /(2M W ), and the loop integration can provide powers of 
the top mass. In the case at hand, the one-loop integrals are convergent, so that the 
one-loop contributions scale at most like g^/M 2 ~ g 2 /M^. At the two-loop level, the 
gaugeless contributions scale as gf/M 2 ~ g 2 g 2 /M^ ~ g 4 M 2 /M^, but the same heavy top 
behaviour can be obtained by inserting a dimension four operator[| proportional to g 2 in a 

1 Dimension two insertions are removed by our choice of renormalization in the W — <j> sector. 
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Figure 3: Insertions of an effective Z° penguin vertex in a one-loop diagram. 

topless loop. In general, the effective lagrangian obtained after integrating out the heavy 



top tells us exactly which the relevant operators are |30|. In our case, only the diagrams in 



Fig. 3 contribute to the leading HTE through the insertion of a flavor changing dimension 
four Z° penguin operator of the kind SLj^b^Z^. The diagram with a mass insertion on 
the internal b line depends on the regularization scheme - it vanishes if IR divergences are 
regulated dimensionally - and in the schemes where it does not vanish it is cancelled in 
the matching by a contribution from the electroweak penguin operator of Eq. In 
both cases, however, its contribution is reintroduced in the quantity C^/ by Cf '. In the 
limit of a heavy top the effective vertex has the form 



with c w = cos9 w . Inserting this gauge-independent effective coupling in the one-loop 
diagrams of Fig. 3, and keeping in mind the tree level couplings of the Z° boson with 6^, 
~ (1/2 — and with bn, ~ — s^/3, we obtain the difference between the HTE of our 

result and the gaugeless limit. The argument is completely analogous for Cg, whose HTE 
also differs from the gaugeless approximation. 

So when does the gaugeless limit potentially fail at two-loop? Whenever at the one-loop 
level the top quark diagrams in the limit of a heavy top scale like a constant, namely in 
the same way as the topless contributions. Indeed, in this case we know that there are 
some dimension four effective operators proportional to gl that can be inserted in one-loop 
diagrams not containing the top and give contributions of the same order, in the limit of 
heavy top, of those belonging to the two-loop gaugeless approximation. 

Table 1 summarizes the situation for the processes considered in the literature at the 
two-loop level. It should be clear that the gaugeless limit works safely only when the 
asymptotic expansion in M t has maximal power (M t 2 at one-loop, M t 4 at two- loop). Of 
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Table 1: Leading HTE contributions to different processes. Following the criterion given in the 
text, the gaugeless limit fails in the first two cases. 



course, there might be exceptions. Indeed, whether the 0{gf) dimension four operators 
are relevant or not depends on the process under consideration. For instance, in the Hgg 
effective vertex [31j - relevant for gluon-gluon fusion production and hadronic decays of 



the SM Higgs boson — they are not because the gluons have no electroweak interaction. 



This is in contrast to the similar case of the if 77 effective vertex [18], where the gaugeless 
approximation does not give the correct result. Similar considerations apply to the heavy 
Higgs limit, although the leading term in the heavy Higgs expansion, subject to other 
constraints, is not always what is expected from dimensional analysis. 

4. Let us now examine in some detail the effect of our calculation on BR 7 . As a first 
step we calculate the Wilson coefficients at a scale fib = 4.8 GeV. It is well-known that the 
large mixing between Q 2 = (cb)v-A {cs)v-A an d Q7 induces additive terms in the running 
of the coefficients from the W to the b mass scale which are numerically very important. 
Our aim is to resum all contributions 0{g 2 a n s L n ) and 0(«M t 2 / 'M^a"L n ), where L is a large 
logarithm. As mentioned above, these terms are uniquely originated by heavy degrees of 
freedom and enter only the determination of the Wilson coefficients at a high scale /1 ~ M w . 
At this order the evolution of the coefficients is therefore driven only by LO QCD effects. 
The Wilson coefficient at the bottom mass scale is given at LO in QCD by 

l(0)eff,,, A _™ r tO)eff, 1ijr \ , 8 / 14 16N (0)e// ^^(0) 



Cr en ^ b ) = v^C^(M w ) + ^ (t/« - Cr"(M w ) + C^(M W ) £ h lV a > (10) 

<J „■ 1 



where 77 = a s (fi w ) / a s (fib) ~ 0.56 and hi, are constants given e.g. in [p9[ . The last term 
is approximately equal to —0.173 and is dominant. 

The electroweak corrections affect Eq. (|T0|) in two ways: (i) they shift Cf )e// (M W ) 



by <JCf{£; (ii) they introduce in Eq. fllUD those gluon and electroweak penguin operators 
that have non-zero 0(g 2 ) or 0(aMf) contributions to the Wilson coefficients at /i = M w . 
In the basis of p9| , these are Q379. Their LO QCD mixing with the magnetic penguin 
operators can be gleaned from the anomalous dimension matrix ^^ e ff given in |i~0|,|34 



after a change of basis. The new entries of ^f^ e ff calculated in [10[] are given in Table 



2 in the conventional basis adopted in ||29|| . The additional contributions to Eq. (|10|) are 
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Table 2: Anomalous dimension matrix entries relevant for the mixing between electroweak and 
magnetic penguins [|1(| in the basis of Two up and three down active flavors are assumed. 



therefore 



5 ew Ci 0)eff (/i b ) =r]%5C e 7 



t=3,7,9 



E 



where Cf (M w ) are the relevant 0(a) contributions to the Wilson coefficients and hij, a.y are 
magic numbers that can be easily determined from the anomalous dimension matrix. The 
last term in Eq. (|TT]) is approximately given by 0.15 C 3 (M W ) + 0.12 C 7 (M W ) - 0.03 C 9 (M W ) 
and is numerically very small; it reduces cf )e// (ji h ) by -0.2%. 

Notice now that Cf\M w ) in Eq. fllCf ) is unaffected by electroweak corrections of the 
kind considered here if is used to normalize the effective Hamiltonian (as in fact we 
do); in that case C2°\m w ) does however receive 0(g 2 s^) corrections. In the NDR scheme 
we find 

w - - • a{M ' 



1 + 



47T 



22 4 , Ml 



(12) 

where a(M w ) is the electromagnetic running coupling evaluated at M w . The 0(a) elec- 
troweak corrections to the Wilson coefficient Ci are therefore very small (—0.13%). But 
we stress that a different choice of normalization would induce additional (much larger) 
electroweak contributions, as can be easily seen using [17]. Eq. ( |P2]) is a new result that 
improves on p5| , where only QED effects were taken into account, and includes all one-loop 
electroweak contributions. 

From Eq. flT0| ) and neglecting the 0(a) effects of Eq. ( |T2"D we see that C 7 ^(fib) is 
reduced by only (1.3 ± 0.2)% for M H = 100 GeV due to the electroweak corrections we 
have calculated; the reduction is less pronounced for larger Higgs masses. In a similar 
way, we find that for Cg (/lb) the numerical impact of the last term in Eq. (O) is more 



important than that of the two-loop correction; electroweak effects increase C^ e ^ (fib) by 
about 0.3 ± 0.2% for M H = 100 GeV. In addition to the leading logarithmic QCD effects 
considered in Eq. ( |10"P there are next-to-leading QCD terms that enter the calculation of 
BR 7 ||. Our electroweak corrections affect them only as a NNLO effect. As there are other 
uncalculated contributions to that order, even in the approximation we have adopted - 
in particular 0(g 2 a s ) corrections to the Wilson coefficients, as can be seen from Eq. (3.14) 
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of [I9| - our corrections should be implemented only in the calculation of the LO QCD 
Wilson coefficient Cy^ "(/it). 

The total effect of electroweak corrections on the NLO calculation of BR 7 is a —2.0 ± 
0.3% reduction for a light Higgs mass, M H w 100 GeV. For larger Higgs masses the effect 
becomes smaller: -1.6% for M H = 200 GeV and -1.3% for M H = 300 GeV. In the 
gaugeless approximation and excluding the last contributions to Eq. (|iT|) , the net effect on 
BR 7 is a 0.5% reduction for M H « 100 GeV. 

We now calculate BR 7 following closely |||,||,||, updating the input parameters and 
introducing only minor refinements in the NLO analysis. For instance, we evaluate the MS 
top mass from the pole top mass using the 0{acf) expression []36[| , as these corrections are 
large and their origin is distinct. Compared to the use of the 0(a s ) conversion formula, this 
leads to a —0.4% reduction of BR 7 . With respect to the detailed analysis of ||, we adopt 
a new CKM factor f CKM = \V t *V tb /V cb \ 2 = 0.97 ± 0.02 instead of 0.95 ± 0.03, obtained 



using < g < 0.4 from global fits of the unitarity triangle [37]. For the semileptonic BR, 
we employ BR5^ = 0.1045 ± 0.0021, corresponding to the T resonance determination (the 
average of LEP measurements is 0.1073 ± 0.0018). We also use M t = M[ ole = (174.3 ±5.1) 
GeV, a s (M z ) = 0.119 ±0.002, M w = 80.419 GeV |!|38|, A 2 = 0.12 GeV 2 , r cb = m c /m b = 
0.29 ± 0.02, M cb = m b — m c = 3.39 ± 0.04 GeV. Employing a conventional definition of 
total BR 7 with a cut on the photon energy E 1 > (1 — S) m b /2, 5 = 0.9, we obtain 

BR - 000399 fcKM BRsL ( M * \^ ( r <* V ™ ( M <*> /iq\ 

BR 7 - 0.000329 — - — J ^— j {—) ^— j . (13) 

Here the dependence of the calculation on the main input parameters is summarized for 
small (< la) variations around their central values. Notice that, compared to the 



2% reduction due to electroweak corrections is compensated by a 2% increase from the 
up-date of the input parameters. The largest present single parametric uncertainty comes 
from r cb and reaches around 5%. 

The choice of 5 = 0.9 for the photon cut-off energy in Eq. ([IB]) is mainly motivated by 
the need to compare with previous literature. The experimental measurement is based on 
a much stronger cut, E 1 > 2.1 GeV [||, but needs to be extrapolated to a more inclusive 
branching fraction (see [f| for a recent discussion). On the other hand, non-perturbative 
problems may arise for soft photons. It seems therefore useful to know BR 7 for higher 
and more realistic cut-offs. For 0.3 < 5 < 0.9 the central value of Eq. ( pT3|) is very well 
approximated by 

BR 7 (5) = 3.01 + 1.0 1 5 - 1.49 5 2 + 0.79 5 3 , (14) 

which shows a mild dependence on the cut-off in a large region of 5 @. 

5. In summary, we have reanalyzed in detail the two-loop electroweak corrections to 
B — > X s 7 and B X s g decays. In order to avoid dealing with presumably small photonic 
effects, in our calculation we have neglected terms proportional to not enhanced by M t 2 . 
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We have also accurately discussed the interplay between electroweak and QCD corrections. 
As a byproduct of the calculation we have presented the complete 0(a) corrections to 
the Wilson coefficient C*2(M W ). The total effect of electroweak corrections on the BR of 
B — > Xs'-f, BR 7 , is a —2.0 ± 0.3% reduction, which is three times larger than in previous 
analyses ]TT|. After inclusion of all NLO QCD contributions, of non-pertubative corrections 



and of all known QED and electroweak effects, we find (for M H = 100 GeV and 5 = 0.9) 

BR 7 = (3.29 ± 0.21 ± 0.21) x 10~ 4 (15) 

where the first error is parametric and based on up-to-date experimental inputs and the 
second one is obtained by scanning the various scales considered in [[| between 1/2 and 
twice their central values and adding 0.3% for unaccounted electroweak higher orders. Had 
we combined the scale ambiguities in quadrature, the second error would have been 0.16. 
For M H = 215 GeV, at the other edge of the preferred region from global fits, the central 
value of BR 7 is 3.30xl0^ 4 . Eq. (|lf|) is in good agreement with the present experimental 
value BR 7 = 3.14 ±0.48 0. 

We are indebted to M. Misiak for pointing out to us a missing term in Eq. (|ll|) and for 
many interesting conversations and to A.J. Buras and P.A. Grassi for helpful discussions. 
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